Gutzwiller density functional calculations of the electronic structure of FeAs-based 
superconductors: Evidence for a three-dimensional Fermi surface 
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The electronic structures of FeAs-compounds strongly depend on the Fe-As bonding, which can 
not be described successfully by the local density approximation (LDA). Treating the multi-orbital 
fluctuations from ab-initio by LDA+Gutzwiller method, we are now able to predict the correct 
Fe-As bond-length, and find that Fe-As bonding-strength is 30% weaker, which will explain the 
observed "soft phonon". The bands are narrowed by a factor of 2, and the d 3z 2_ r 2 orbital is pushed 
up to cross the Fermi level, forming 3-dimensional Fermi surfaces, which suppress the anisotropy 
and the (n, n) nesting. The inter-orbital Hund's coupling J rather than U plays crucial roles to 
obtain these results. 

PACS numbers: 74.25.Jb, 71.27.+a 



The iron pnictides are interesting not only because 
of the high superconducting transition temperature Tc 
(above 50K 1]), but also because of their different as- 
pects compared with the high Tc cuprates: mainly (1) 
iron pnictides are typical multi-orbital systems, where 
spin, orbital and charge degrees of freedom are all ac- 
tive 0, while in cuprates an effective single band model 
can be established; (2) the electron correlation strength 
is intermediate 0, 0,0 > n °t as strong as that in cuprates. 
However as will be shown in this paper, due to the multi- 
orbital nature, where the inter-orbital interaction be- 
comes important, the correlation in iron pnictides still 
plays crucial roles in determining not only the correct in- 
ternal structure but also the correct electronic structure 
near the Fermi Surface (FS). 

Both the 1111 (such as LaOFeAs 0) and 122 (like 
BaFe2As2 p|) systems contain the FeAs- layers as the 
building blocks, and the low energy bands around the 
FS are mostly from the five Fe-3d orbitals @, El, 0, 0] 
hybridized with the As-4p orbitals. Since the correla- 
tion is not so strong, the LDA level calculations success- 
fully predict two basic aspects: (1) the metallic electronic 
states 0]; (2) the FS nesting and the existence of SDW 
state in parent compounds @, However, LDA cal- 
culations encounter problems with more and more ex- 
perimental information accumulated. The most serious 
problem is that LDA can not describe the Fe-As bond- 
ing accurately (with error bar as large as 10%) fldi l, but 
the electronic structures of these compounds depend on 
the Fe-As bonding sensitively fill ], leaving the calculated 
electronic structures so far questionable. Since the possi- 
ble pairing mechanisms strongly depend on the detailed 
electronic structures, the correct description to the Fe- 
As bonding and hence the electronic structure is highly 
desirable, which is the main subject of our present paper. 

The fact that LDA band-width is about two times 
wider than ARPES measurement [H, 0, 0] indicates that 



the fail of LDA is mostly due to the insufficient treat- 
ment of electron correlation. Unfortunately, the fully self- 
consistent treatment of multi-orbital fluctuation is still a 
challenging task up to now. The correlation effects in 
iron pnictides have been studied by several groups [l2j 
using DMFT (dynamical mean field theory). Although 
the structure optimization is not achieved, these studies 
indicate that they are in the intermediate coupling re- 
gion with the effective mass enhancement varying from 
1.5 to 3.0. Using our newly developed LDA+ Gutzwiller 
method [131 ] , where density functional theory is combined 
with the Gutzwiller variational approach with full charge 
self-consistency, we are able to do structure optimizations 
with orbital fluctuations included. It has been shown [l3| 
that, for the ground state determination, the energy ac- 
curacy of Gutzwiller approach is comparable to DMFT. 
We will show in this paper that Fe-As bonding can be well 
described after considering the correlation: (1) the Fe- 
As bond length can be correctly predicted; (2) the Fe-As 
bonding strength is now estimated to be 30% weaker than 
that of LDA, which explains the experimentally observed 
soft phonon modes [3| • We realized that to get those cor- 
rect understanding, the Hund's coupling J rather than 
Coulomb U is the crucial factor for such multi-orbital 
systems. 

With the correct Fe-As bonding obtained, we study 
the normal state electronic structure, and show that FS 
are quite different with previous LDA results. The ef- 
fect of correlation will shift the d 3z 2_ r 2 orbital up to the 
Fermi level, as the result 3-dimensional (3D) FS appear 
in the normal state. The existence of 3D FS was im- 
plied by many experiments, i.e. the low anisotropy in 
normal state resistivity, up critical field and superfluity 
density 15, 16|, 17, 18| . Moreover, from the calculations 
of Linhard function, we found that such careful treatment 
of orbital fluctuation will reduce the (ir, n) FS nesting ef- 
fect as well as the effective moment in SDW phase. 
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FIG. 1: The optimized (by LDA+G) Fe-As inter-layer dis- 
tance dp e _ As and band-width renormalization factor Z 2 as 
function of interactions (U and J). The left and the right 
panels are the results for LaOFeAs and BaFe2As2 respectively. 



The calculations are done with the BSTATE (Beijing 
Simulation Tool for Atom Technology) code, in the ultra- 
soft pseudopotential plane wave method (l3j . To guaran- 
tee the convergence, we use 30 Ry for the cut of energy 
of wave-function expansion, and 12x12x8 K-points mesh 
in the Brillouin zone. The projected wannier functions 
are used as the local orbitals, and the virtual crystal ap- 
proximation (VCA) is used for doping. Since our local 
basis is wannier function with all states included, it is 
therefore a reasonable value of £/=3.0-4.0eV, and J=0.8- 
l.OeV, as suggested by Anisimov Nevertheless, differ- 
ent parameters has been applied, and optimized values 
(Z7=3.0eV, J=0.8eV) can be obtained from the correct 
structure determination. 

1. Fe-As bond length, phonon frequency, and band nar- 
rowing. 

The equilibrium position of As, described by d z Fe _ As 
(Fe-As inter-layer distance along the c axis) is studied, 
and the results are shown in Fig. 1(a) and (b). The LDA 
d z Fe _ As is 0.15A (for LaOFeAs) and 0.16A (for BaFe 2 As 2 ) 
shorter than their corresponding experimental value, this 
is in agreement with previous calculations [Hi ]. However, 
with the increasing interaction strength, the d Fe _ As in- 
crease, and approach the experimental value with the pa- 
rameters U ^3.0eV and J ~0.8eV. The same parameters 
reproduce the As equilibrium position for both LaOFeAs 

is sensitive to J rather 



and BcFe2As2, and the d 



Fe-As 

than U, implying the important role of J. The same 
calculations have been repeated for LaOo.gFo.iFeAs and 
Bao.6Ko.4Fe2As2, and the same results are obtained. 

Not only the Fe-As bond-length but also the Fe-As 
bonding strength are seriously affected by the correla- 
tion. It has been a puzzling issue that the measured 
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FIG. 2: (a) The total energy versus As-displacement (Ad) 
from its equilibrium position, (b) The occupation numbers of 
five Fe-3d orbitals as function of J. U is fixed to be 3.0eV. 
(c) The projected density of states (PDOS) of Fe-3d 3z 2_ r 2 
orbital for different situations: LDA for bulk, LDA+G for 
bulk, and LDA for surface. Please note the relative change of 
3d 3z 2_ r 2 weight at the Fermi level, (d) The Lindhard response 
functionn Xo(<?) f° r g- vector along the F — M line. 



phonon spectrum of iron pnictides show certain soften- 
ing compared with LDA results [l4|. There are three 
main peaks in the phonon density of states (DOS), lo- 
cated around 13, 23, 30 meV, respectively [14[. The two 
peaks (at 13 and 23meV) can be reproduced from LDA, 
however the peak around 30mcV, which are mostly due 
to Fe-As bond related modes, is 15% softer than LDA 
results. This "soft phonon" energy can be obtained after 
reducing the Fc-As bonding strength (force constants) by 
30% from its LDA values artificially. Here we show that 
this weakening is naturally explained by the correlation 
effect. Fig. 2(a) shows the calculated total energy as func- 
tion of As-displacement. The curves can be nicely fitted 
with parabolic equations, suggesting the minor anhar- 
monical effect. From the second derivative of the curves, 
we can estimate the force constant, which is 30% smaller 
in LDA+G comparing with LDA. This reduced force con- 
stant will exactly reproduce the observed "soft phonon" 
around 30meV. We have also optimized the La-position 
in the same way, and found that its equilibrium posi- 
tion and force constant are little affected by the LDA+G 
treatment. 

It has been argued in literature that the discrepancies 
discussed above can be improved if spin-polarization is 
assumed [20J, while our strategy is to treat the system 
without such assumption. This is of course the true ex- 
perimental situation, where the normal state above tran- 
sition temperature is paramagnetic metal. In the pres- 
ence of orbital fluctuation, the essence of Gutzwiller ap- 
proach is to suppress the weight of double-occupancy, 
and renormalize the kinetic energy by the factor Z 2 (av- 
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FIG. 3: The electronic structure of LaOo.gFo.iFeAs. Upper 
and middle panels: Band Structures by LDA and LDA+G 
respectively. We plot the fat bands, where the symbol size 
corresponds to the projected weight of Bloch states onto the 
Fe-3d 3z 2_ r 2 orbital. Lower panels: Fermi surfaces from LDA 
(left) and LDA+G (right). The F point located at the body 
center of the cubes. 

eraged through five orbitals) as shown in Fig. 1(c) and 
(d). For the U and J determined from structure opti- 
mization, the renormalization is about 0.6, which is con- 
sistent with the ARPES results 0, 0, 0] . (The Z 2 factor 
is again not sensitive to U, but sensitive to J). The nar- 
rowed band-width will favor the weaker Fe-As bonding 
and longer bond-length. It is worth to note that assuming 
spin polarization is in fact the simplest way to suppress 
the double-occupancy although it is no longer realistic 
for the normal state of superconducting compounds. 
2. The role ofd 3z 2_ r 2 orbital. 

It is essential to understand the role of Hund's cou- 
pling J and discuss the change of orbital characters in 
the presence of interactions. Starting from the struc- 
ture aspect, the crystal splitting of d orbitals is not so 
strong (of the order 0.1-0.2eV), nevertheless the squeez- 
ing of FeAs4 tetrahedral along c-axis push the d 3z 2_ r 2 
orbital down to the lowest energy among five, and it con- 
tribute little to the DOS at Fermi level (Ef). In most 
of the effective models studied so far, this orbital is ne- 
glected for simplicity. In the presence of interactions, 
the role of Coulomb U is to enhance the orbital polariza- 
tion, however the role of inter-orbital Hund's coupling J 
is opposite: it try to tight all orbitals together and favors 
even distribution of the 6 electrons among five orbitals. 
As shown in Fig. 2(b) and (c), with the increasing of J, 
the energy level of d 3z 2_ r 2 is raised and its occupation 
is reduced. For the optimal J=0.8eV, all orbitals are al- 
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FIG. 4: The electronic structure of Bao.6Ko.4Fe2As2. All no- 
tations are the same as Fig. 3. For the comparison reason, 
we plot the band structures in the same high symmetry lines 
as that in LaOFeAs, instead of the BZ of the body centered 
structure. 

most equally occupied. The reduced orbital polarization 
will enhance the inter-orbital fluctuation for such multi- 
orbital system. This is the key to understand the kinetic 
energy renormalization factor Z 2 around 0.6. 

The d 3z 2_ r 2 orbital is now relevant to the low energy 
physics, very different with previous understanding based 
on LDA. Most significantly, 3D FS appears and the (%, ir) 
FS nesting is suppressed. The calculated Linhard re- 
sponse funciton shown in Fig. 2(d) suggest that the nest- 
ing is not as strong as in LDA, and slightly incommensu- 
rate. The SDW solution with smaller magnetic moment 
(comparing with LDA) would thus be expected for the 
parent compounds. 

3. The band structures and FS. 

The normal state band structures of superconducting 
compounds, LaOo.gFo.iFeAs and Bao.6Ko.4Fe2As2, are 
shown in Fig. 3 and 4. The calculations are all done with 
the experimental As-position, rather than the incorrect 
position optimized from LDA. In the LDA results, the 
bands which has strong d 3z 2_ r 2 component are away from 
the Ef, and the FS are highly 2-dimensional. In the 
LDA+G results, however, those bands are clearly shifted 
up and go cross the Fermi level. As the results, one of the 
FS is strongly 3D due to the dominate d 3z 2_ r 2 character. 
The enhanced d 3z 2_ r 2 contribution to the DOS at Ef can 
be quantitatively obtained from Table. I. 

4-. Comparison with experiments. 

Although most of the iron pnictides have layered struc- 
ture, the anisotropy reflected in transport and super- 
conducting properties is quite low. For instance, in 
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TABLE I: Some calculated parameters: ^f^- is the conduc- 
tivity ratio between a and c, N 3z 2_ r 2 is the partial DOS of 
d 3z 2_ r 2 orbital, and N tota i is the total DOS, at the Fermi 
level. 
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Ba(Fen.32HCon. 074)2 As2, the anisotropy in resistivity is 
around 3±1 [15j. the anisotropy in penetration depth is 
around 3 to 6 [lg, [l7|, and that of the up critical field 
for hole doped compound Bai_ x K x Fe2As2 is around 2 
to 3 just below the critical temperature but down to 
unity around 10K |18l |. The inelastic neutron scattering 
measurement also show strong 3D features in the reso- 
nance peaks appearing in the superconducting phase [2l| . 
Those low anisotropy results can be naturally explained 
by the existence of the 3D FS, which appears only when 
the correlation effect has been carefully included with 
the correct experimental structure. We have calculated 
the resistivity anisotropy based on the band structure 
obtained by LDA+G assuming the isotropic relaxation 
time. The results summarized in the Table I fit quite 
well with the experimental results. In addition, the pen- 
etration depth [3|j_5Yt 0' an d thermal conductiv- 
ity measurements [24] imply there may be line nodes in 
the superconducting states. From our results, we pro- 
pose that, at least for the Bao.6Ko.4Fe2As2, the 3D FS at 
k z — 7r plane is big enough to cross the node lines of the 
S± pairing state [9(. Therefore the line nodes may exist 
on this hourglass like 3D FS in the S± pairing state. 

On the other hand, however, the ARPES results @, 
suggest discrepancies with transport measurements: (1) 
all the detected four FS are cylinder like with very weak 
dispersion along the c-axis; (2) the superconducting gaps 
in all the four FS are isotropic without any line nodes. 
It is still too early to address these discrepancies, while 
one possible scenario based on our calculations is the fol- 
lowing. Since the main component of the hourglass like 
3D FS is d^ z 2_ r 2 1 which disperses strongly along the c- 
axis, it should be strongly affected by the surface. To see 
the general tendency of surface effect, we have done the 
surface calculation using LDA only. We can see from 
Fig. 2(c), the PDOS of d 3z 2_ r 2 at Ef is strongly sup- 
pressed on the surface, which makes it very difficult to 
be detected directly by ARPES. Nevertheless, the possi- 
ble k z dependence of FS has been partly observed in the 
recent ARPES measurement [25j |. 

In summary, by properly treating the electron correla- 
tion through self-consistent LDA+Gutzwiller method, we 
are now able to describe the Fe-As bonding (both bond- 
length and bonding-strength) successfully. Having this 
prerequisite, the resulting electronic structure is very dif- 



ferent with previous LDA understanding. In particular, 
the bands are narrowed by a factor of 2, and the d 3z 2_ r 2 
orbital is now relevant to the low energy physics by con- 
tributing to a 3D FS, which suppress the anisotropy and 
the (it, 7t) FS nesting. Although the interaction strength 
in iron pnictides is not as strong as in cuprates, the inter- 
orbital Hund's coupling J (due to the multi-orbital na- 
ture) plays crucial roles to determine the electronic struc- 
ture. 
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